Dense bubble flow in a silo: an unusual flow of a dispersed medium 
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The dense flow of air bubbles in a two-dimensional silo (through an aperture D) filled with a 
liquid is studied experimentally. A particle tracking technique has been used to bring out the main 
properties of the flow: displacements of the bubbles, transverse and axial velocities. The behavior 
of the air bubbles is observed to present similarities with non-deformable solid grains in a granular 
flow. Nevertheless, a correlation between the bubble velocities and their deformations has been 
evidenced. Moreover, a new discharge law (Beverloo-like) must be considered for such a system, 
where the flow rate is observed to vary as D 1 ^ 2 and depends on the deformability of the particles. 

PACS numbers: 47.50.+d, 47.60.+i, 82.70.-y, 83.50.-v, 83.80.-k 



Structured fluids like granular materials or aqueous 
foams received much attention from physicists during the 
last decade. Indeed, they are ubiquitous in nature and 
in industrial processes Q, Q, and exhibit striking The- 
ological behaviors as compared to conventional liquids. 
These materials involve multiple physical processes and 
cooperative phenomena such as shear banding, particle 
interactions and the formation of arches redirecting the 
forces inside the material 0, • Such processes are then 
able to provoke intermittent flows Q or jamming Q be- 
cause of the topological constraints that develop since 
the particles (grains or bubbles) press against each other. 
This contributes to the complex flow behavior of granular 
media and foams. 

The present study deals with bubble assemblies where 
the bubbles are predominantly oblate so that the contacts 
between two adjacent bubbles are reduced (see Fig.njj). 
The liquid fraction in the material is high enough to con- 
sider the diffusion of air through the bubble boundaries 
as negligible at the time scale of the experiments. 

Bubble assemblies can also be seen as a granular 
medium where contacting particles are characterized by a 
nearly zero friction. Thus, we have a remarkable system 
with a low energy dissipation at the contacts in com- 
parison with traditional grain assemblies. Both granular 
media and foams are model systems for divided materi- 
als. Nevertheless, even if bubble flows present qualitative 
analogies with granular flows, we must keep in mind that 
the nature of the interactions between the particles is to- 
tally different. Moreover, in contrast with granular mate- 
rials, the bubbles may undergo strong deformations when 
they are subjected to a constraint. Many works deal with 
the study of flowing foams (velocity profiles, bubble de- 
formations) in a Couette geometry 0, H or through a 
constriction 9J. In contrast with this latter paper deal- 
ing with a dry foam flowing in a Hele-Shaw cell, our in- 
terest concerns wet foams flowing between a plate and 
the liquid surface. Instantaneous measurements of veloc- 
ity and deformation of the bubbles are reported. The 
dependence of the bubble flow rate regarding the aper- 
ture or the mean bubble velocity is analyzed through the 
comparison with granular materials in a same geometry. 



In this paper, the flow of a two-dimensional (2D) foam 
is analyzed experimentally 0, [n| . The setup consists of 
a transparent inclined glass plane of width W = 130 mm 
which is immersed into water (Fig. QJi) . The tilt angle 8 
of the plane can be adjusted by a fine screwing system. 
Quasi-monodisperse air bubbles are then injected at the 
bottom of the inclined plane by blowing air through a 
needle (the height h of the piling is kept roughly constant 
during the experiments h ~ 200 mm). Note that the tilt 
angle is small (0 ranges from 0° to 1°) to ensure that 
only one layer of bubbles is created in the perpendicular 
direction to the plane and to decrease the influence of the 
gravity (i.e. liquid drainage). Bubble size is controlled by 
an air pump and is kept roughly constant in the present 
study [diameter d = (5.4±0.3) mm]. In order to avoid the 
coalescence of the bubbles, dish- washing liquid based on 
anionic surfactant AEOS-2EO is mixed into water [l2T |. 
The surface tension of the solution is 7 ~ 25mNm _1 
leading to a capillary length k^ 1 = y^f/ifig) — 1.6mm. 
In our experiment, the bubble are larger than kT x and 
consequently adopt an oblate shape (Fig. QJd) so that the 
contact between two adjacent bubbles is reduced. Due to 
buoyancy, the bubbles rise underneath the inclined plane 
and tend to pack on a transverse wall placed at the top 



inclined 
plane 




FIG. f : (a) Sketch of the experimental setup. A transparent 
inclined plane is immersed into water and tilted at an angle 
8. Small air bubbles are injected from below at the bottom of 
the plane and rise towards the top where they aggregate. Mo- 
tions of the bubbles are recorded by a CCD camera, (b) Side 
view: the bubbles present an oblate shape, (c) Top view: the 
bubbles have a circular shape in the bulk. 
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FIG. 2: Transverse v x and longitudinal v y velocity profiles 
during the discharge of a silo (9 ~ 0.63°, D ~ 8.6 mm). (- 
-) position (xq, yo) of the center of the outlet of the silo. 
The different curves (▼, •, ■. A) correspond respectively to 
velocity profiles at increasing distances from yo = (Figs. 
and ^jp) and increasing distances from xo (Figs. and|2Ji). 



of the plane. An orifice of width D at the center (xq, yo) 
of this obstacle allows the bubbles to empty out the silo. 
Our experiment is then analogous to the gravity-driven 
granular flow in a flat bottomed 2D silo, emptying out of 
a central orifice ^| . Top views of the bubble packing are 
recorded through the transparent tilted plane by means 
of a CCD camera at a frame rate of 50 frames per second, 
fn order to quantify and to measure the main properties 
of the flow, the motion of each bubble (ss 2000 bubbles 
per frame) has been tracked through image analysis. 

Figure |3 displays typical transverse v x and longitudi- 
nal v v velocities of the bubbles as a function of x and 
y at different heights in the packing, fn the top part 
of the silo [(A) in Fig. [2b], a block motion is observed: 
the longitudinal velocity v v remains constant reflecting 
low interactions between the bubbles and the side walls 
(slipping motion). As y decreases {i.e. the bubbles move 
towards the outlet), the flow takes place in a triangular- 
shaped domain where the fastest bubbles are located in 
the center of the silo and the slowest or motionless ones 
at the walls [(T) in Fig. [2b]- This is in agreement with 
the 'V-shape' of mobile grains observed during the flow 
of granular materials in silos or hourglasses jl3L Il4| . Far 




FIG. 3: (Color online) Typical recordings of the bubble dis- 
charge: (a) Snapshot of the experimental flow, (b) Displace- 
ment held. The arrows point out a dislocation and a recircu- 
lation zone, (c) Velocity field. The fastest bubbles appear in 
blue (or darker) . (d) Deformation held. The more constrained 
bubbles appear in red (or darker) while circular bubbles ap- 
pear in blue. 



from the outlet (y > 15d), bubbles have a nearly zero 
transverse velocity (Fig. and a constant longitudinal 
velocity v y (Fig. |5Ji). For lower values of y, the velocity 
distribution of the bubbles deep inside the silo is totally 
modified: a transverse component of velocity v x appears 
in the flow conducting the bubbles towards the orifice 
while v y increases strongly near the outlet. 

A zoom of the bottom part of the silo of height 60 mm is 
shown in Fig. [3k.. Image analysis allows one to extract the 
position of each bubble in the stack and evaluate the dis- 
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placement profile (Fig. [2)) and the velocity field (Fig. Ot) 
during the flow. Figure 01 displays the deformation field 
and provides information concerning the constraints un- 
dergone inside the dense bubble assembly. In contrast 
with granular materials where the flow occurs only in a 
cone-shaped central region of the silo (while the grains 
located in the regions near the side walls are motionless) 
13 1 note that motions of bubbles are detected ev- 
erywhere in the silo. Moreover, the bubble tracking put 
into evidence recirculation zones and the propagation of 
dislocations (shear bands) during the bubble flow. As 
expected in two-dimensional nearly monodisperse flows, 
bubbles tend to pack in an ordered hexagonal structure at 
the top of the bubble packing (Fig. |2^) ; these organized 
domains are separated by the dislocations. Motions of 
blocks of bubbles are observed and look like crystal do- 
mains. Those blocks are sliced by fast-moving defects 
along bubble lines. The movement is not simultaneous 
along the whole row but begins at one end with the ap- 
pearance of a dislocation running along the slip line. This 
process is initiated at the orifice of the silo and propa- 
gates towards the bulk after many rebounds on the side 
walls. Both recirculation zone and dislocation appear 
clearly as pointed out by the arrows in Fig. ^p. 

The instantaneous velocity v = (v%, + v 2 ) 1 / 2 of each 
bubble has been computed and superimposed on the im- 
ages of the flow (the darkest - or blue - bubbles corre- 
sponding to the fastest ones). The recirculation regions 
and the dislocation observed in Fig.|3jD are characterized 
by a mean velocity higher than in the rest of the silo 
(Fig. HO- 

The deformation of a bubble is evaluated by fitting the 
bubble shape by an ellipse and computing its eccentricity 
[15). As pointed out in Fig. [3J1, the fastest bubbles are 
observed to correspond to the less deformed ones. Mov- 
ing bubbles press against surrounding ones and deform 
them. This is especially the case near the dislocation. 
Therefore, in such dense bubble flows, spherical bubbles 
are located far from the outlet and propagate in a block 
motion at a constant velocity or in the dislocation, while 
flattest ones correspond to crushed bubbles against the 
bottom wall or constraint bubbles inside the pile. 

A classic paper on the flow of particles through orifices 
is that of Beverloo et al. ^(|. Using a straightforward 
dimensional analysis, they pointed out that the grain flow 
rate Q g is a power law in D (where D is the diameter of 
the orifice). The power depends on the dimensionality of 
the flow [IlH[l7|, so that for a 2D flow Q g oc g 1 / 2 D 3 / 2 . 
A corrective term is usually added to take into account 
the mean diameter of the grains d g : 

Q g <xg l l 2 {D-kd g f/ 2 , (1) 

where k is an empirical constant depending on the parti- 
cle shape and the grain interactions (fc ~ 1.5 for beads). 
This law characterizing the dynamical properties of the 
grain flow does not invoke explicitly dissipative mecha- 
nisms. However, we expect that k and the value of the 




5 10 15 20 25x10-' 



FIG. 4: Bubble flow rate Q as a function of the tilt angle 6 of 
the plane. The curves correspond to different apertures D of 
the silo: (T) D = 6.3 mm, (■) D = 8.0 mm, (•) D = 9.7mm, 
(♦) D = 12.8 mm, (a) D = 19.9 mm; (— ) best fit of the 
experimental data with Q oc (sin^) 3 '' 2 . 



exponent 3/2 contain implicitly friction and deformation 
of the grains even if the dependence of these parameters 
with friction and deformability is still an open question. 

The bubble flow rate Q has been evaluated by counting 
the number of bubbles flowing out of the silo, at different 
tilt angles 8 of the plane. Figure 0] shows that Q varies 
as (sin6*) 3 / 2 , where the tilt angle plays a similar part 
to gravity encountered in the traditional Beverloo's law. 
The behavior of bubble flow through an aperture is thus 
drastically different from the one of a granular flow for 
which Q g varies as g 1 ^ 2 - This is probably due to the 
nature of forces between the particles and with the walls: 
the frictions inside a granular medium are much more 
important than in a bubble flow and consequently slow 
down the silo discharge. 

Figure [3] displays the bubble flow rate Q as a func- 
tion of the normalized aperture of the silo D/d. In con- 
trast with granular materials for which the flow rate Q g 
vanishes for D/d < 1.5 due to the formation of arches 
[Eq. ||TJ], Q is observed to be non-zero for smaller values 
of the aperture. Bubbles can pass through an orifice of 
the order of 2<i/3 after having undergone strong shape 
deformations. The value of k ~ 2/3 should then de- 
pend mainly on the bubble size and the surface tension 
of the liquid; these assumptions will receive special care 
in a subsequent work. Moreover, note that Q varies as 
(D/d) 1 / 2 which is significantly different from the case of 
a non-deformable 2D granular material flowing down a 
silo where the power is 3/2. 

The bubbles at the outlet of the silo are affected by 
the buoyancy B coming from the bubbles located behind 
them. This buoyancy presents a non-trivial law because 
it does not result directly from the sum of each individual 
bubble in the stack: a part of the energy is dissipated 
through the deformation of the bubbles. Let assume that 
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FIG. 5: Bubble flow rate Q as a function of the aperture 
D of the silo (normalized by the mean bead diameter d). 
The curves correspond to different tilt angles 9 of the plane: 
(T) 9 = 0.46°, (o) 9 = 0.50°, (■) 9 = 0.56°, (o) 9 = 0.63°, 
(•) 9 = 0.75°, (A) 9 = 0.82°; (— ) best fit of the experimental 
data with Q oc (D/d — k) 1 ^ 2 ; k presents very low dispersion 
and verifies k = (0.66 ± 0.03) ~ §. 



the driving force is the buoyancy component along the 
plane, so that B oc g&inO. 

The dissipation relationship linking the drag force / 
to a given velocity v is non-tri vial, and known to involve 
rather subtle hydrodynamics 0, [l^. If Archimedes' 
forces were balanced by a Stokes' viscous force (/ oc r]dv ), 
this would provide a constant rise velocity proportional 
to sin 6. But for a single bubble creeping along a slightly 
inclined plane immersed in a viscous fluid, it has been 
shown that the creeping velocity v does not vary linearly 
with the tilt angle sin 9 [2(j. Moreover, the viscous dis- 
sipation is known to be sensitive to the surface mobility 
of the surfactant For our surfactant which leads to 
tangentially mobile bubble surfaces, if we consider the 
viscous dissipation associated with the formation of a lu- 
bricating film between the single bubble and the plane, 
the viscous force / is observed to vary as 7 1 / 12 (7yi;) 2 / 3 . 

Assuming in a first approximation that the dynamics 
of the flow is mainly governed by these two forces (ne- 
glecting therefore the interactions between the bubbles), 
the equilibrium between the buoyancy and the drag force 
leads to a terminal velocity v proportional to (sin#) 3//2 . 
Therefore, the bubble flow rate Q = vD/d 2 can be writ- 
ten 

(g sin 0) 3 / 2 
Qrx V1 i/s ' ( 2 ) 

which agrees with behaviors emphasized in Fig. 01 More- 
over, as expected, the bubble flow rate is then observed 
to decrease when the viscosity r\ of the fluid is increased. 
This straightforward model reproduces the main effects 
observed experimentally and enhances the fact that a 
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FIG. 6: Rescaling of all measurements using Eq. © for the 
data presented in Figs. 0] and 03 



new Beverloo's law must be considered in the particu- 
lar case of bubble flows and more generally the flow of 
deformable dispersed media. 

Figure [fj] displays a rescaling of all measurements and 
confirms that the Beverloo's law for a deformable dis- 
persed medium is given by: 

Q o, (g sin9 f 2 (^-k^j ' , (3) 

where fc < I is related to the deformation ability of the 
particles and the friction between particles. More exper- 
iments must now be performed to distinguish the influ- 
ence of both friction and deformation on k 7 and on the 
value of the exponent. Numerical simulations are cur- 
rently performed on this topics and will be the object of 
a subsequent publication. Finally, following the approach 
proposed by Asipauskas et al. |9j , the study of the influ- 
ence of the bubble deformation on the flux measured at 
the outlet of the silo would provide a good complement 
to the present study. 

In summary, the present experiment shows that the ve- 
locity profiles for bubble flows in silos present qualitative 
similarities with a granular flow but with sliding motions 
at the walls. A correlation between the velocity and the 
deformation of the bubbles has been observed. More- 
over, we demonstrate that a new Beverloo's law must be 
considered for the specific case of deformable dispersed 
media flowing out a silo, where the flow rate varies as 
D 1 ! 2 and depends on the nature of the interactions be- 
tween particles and their capability to deform themselves 
when they are submitted to a constraint. 

The authors thank G. Broze and S. Dorbolo for helpful 
discussions during this work. Y. B. benefits from a post- 
doctoral fellowship supported by the ARC programme 
number 02/07-293 and ESA MAP "AO-99-108". 
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